Vortex clusters and multiquanta flux lattices in thin films of 

anisotropic superconductors 

A. V. Samokhvalov«, D. A. Savinov^, A. S. Mel'niW 1 ) , A. I. Buzdin^ 

W Institute for Physics of Micro structures, 

Russian Academy of Sciences, 
603950 Nizhny Novgorod, GSP-105, Russia 
( 2 ) Institut Universitaire de France and Universite Bordeaux I, France 

Abstract 

The distinctive features of equilibrium vortex structures in thin films of anisotropic supercon- 
ductors in tilted magnetic fields are studied for the limits of moderate and strong anisotropy. The 
energetically favorable shape of isolated vortex lines is found in the framework of two particu- 
lar models describing these limiting cases: London theory with an anisotropic mass tensor and 
London-type model for a stack of Josephson-decoupled superconducting layers. The increase of 
the field tilting is shown to result in qualitative changes in the vortex-vortex interaction poten- 
tial: the balance between long-range attractive and repulsive forces occurs to be responsible for a 
formation of a minimum of the interaction potential vs the intervortex distance. This minimum 
appears to exist only for a certain restricted range of the vortex tilting angles which shrinks with 
the decrease of the system anisotropy parameter. Tilted vortices with such unusual interaction 
potential form clusters with the size depending on the field tilting angle and film thickness or/and 
can arrange into multiquanta flux lattice. The magnetic flux through the unit cells of the corre- 
sponding flux line lattices equals to an integer number M of flux quanta. Thus, the increase in 
the field tilting should be accompanied by the series of the phase transitions between the vortex 
lattices with different M. 

PACS numbers: 
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I. INTRODUCTION 



According to a standard picture of the mixed state in bulk type-II superconductors the 
Abrikosov vortices penetrating the homogeneous sample form a periodic arrangement called 
a flux lattice 1 . The magnetic flux through the unit cell of such flux line lattice equals to the 
flux quantum 0o — nhc/e: we have one vortex per unit cell. There are a few examples of 
rather exotic superconducting systems which may provide a possibility to observe a differ- 
ent vortex lattice periodicity, namely the structures with more than one vortices per unit 
cell. In particular, the phase transitions to such multiquanta flux lattices can occur, e.g., 
for superconductors with unconventional pairing^ or 2D Fulde-Ferrell-Larkin-Ovchinnikov 
superconductors*. 

The goal of this work is to suggest an alternative scenario of the phase transitions between 
the flux structures with different number of vortices per unit cell which can be realized 
in thin films of anisotropic superconductors. The underlying physical mechanism for this 
scenario arises from the interplay between the long range attraction and repulsion between 
tilted vortex lines in thin films discussed recently in Ref. Ill The unusual attractive part of 
the vortex-vortex interaction potential is known to be a distinctive feature of anisotropic 
superconductors and the value of the attractive force is controlled by the tilting angle of the 
vortex line with respect to the anisotropy axis 6-8 . The origin of the long range intervortex 
repulsion in thin films has been analyzed in the pioneering work^ by Pearl in 1964. This 
repulsion force always overcomes the attraction at rather large distances because of the 
different power decay laws of these contributions. Note that, of course, the short range 
interaction between vortices is also repulsive. Finally, this balance between the repulsion 
and attraction can result in the formation of the nonmonotonic interaction potential U(R) 
vs the intervortex distance R. Increasing the vortex tilting angle we first strengthen the 
attraction force between vortices and, thus, the minimum in the vortex interaction potential 
can appear only for rather large tilting angles when the attraction overcomes the Pearl's 
repulsion. This minimum shifts towards the larger intervortex distances with the further 
increase in the tilting angle and, finally, at rather large distances the attraction appears to 
be suppressed due to the exponential screening effect. As a consequence, the minimum in 
the interaction potential exists only for a certain restricted range of the vortex tilting angles 
which shrinks with the decrease of the system anisotropy parameter. The appearance of a 
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minimum in the interaction potential points to the possibility to get a bound vortex pair (or 
even the clusters with higher vorticities) for a certain range of vortex tilting angles. For a 
flux line lattice such vortex-vortex interaction potential can cause an instability with respect 
to the unit cell doubling, i.e. the phase transition to the multiquanta vortex lattices. 

In this paper we use two theoretical approaches to describe the peculiarities of the inter- 
vortex interaction and resulting formation of clusters and multiquanta lattices. One of them 
is a standard London model accounting for an anisotropic mass tensor which is adequate 
for the superconductors with moderate anisotropy. This approach assumes that the super- 
conducting coherence length in all directions exceeds the distance between the atomic layers 
and obviously breaks down in the limit of strong anisotropy, i.e., for Josephson-coupled 
layered structures. In the latter case we choose to apply another phenomenological model, 
namely the so-called Lowrence-Doniach theory™. For rather small intervortex distances this 
theory can be simplified neglecting the effects of weak interlayer Josephson coupling. This 
approach of Josephson-decoupled superconducting layers is known to be useful in studies of 
the vortex-lattice structure at low fields^ 1 ^. 

Considering thin film samples in tilted magnetic fields we do not restrict ourselves by 
the case of only straight vortex lines and study the problem of the energetically favorable 
vortex line shape in the presence of the inhomogeneous supercurrent screening the field 
com pon ent Hn parallel to the film plane. Previously this problem has been addressed in 
for rather small deviations of the vortex line from the direction normal to the film 



Ref. 
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plane. Such approximation is obviously valid only for the | Hn | values much smaller than 
the critical field of the penetration of vortices parallel to the film plane. For anisotropic 



London model this analysis of Ref. \l4 h as been previously generalized for the case of a 



strongly distorted vortex line (see Ref. Il5l). For the sake of completeness we present here the 



calculations of the shape of an isolated vortex line for arbitrary fields | Hi| < within 
both theoretical models describing the limits of strong and moderate anisotropy. As a next 
step, we calculate the vortex-vortex interaction potential for such strongly deformed vortex 
lines. Further analysis in the paper includes the calculations of energy of finite size vortex 
clusters as well as the energy of vortex lattices with different number of vortices per unit 
cell. 

Experimentally the visualization of unconventional vortex arrangements could be car- 
ried out by a number of methods which provided convincing evidence for the existence 
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of vortex chains in bulk anisotropic superconductors caused by the intervortex attrac- 
tion phenomenon (such as the decoration technique in Y BaiCuj,0'r^- 1 scanning-tunneling 
microscopy in NbSe^, scanning Hall-probe^ and Lorentz microscopy measurements in 



The paper is organized as follows. In Sec. II we find the energetically favorable shape of 
an isolated vortex line. In Sec. Ill we calculate the vortex-vortex interaction potential and 
prove the existence of a potential minimum for a certain range of field tilting angles and 
parameters. The Sec. IV is devoted to the calculation of energy of vortex clusters. Finally, 
in Sec. V we present our analysis of the phase transition between the vortex lattices with 
one and two flux quanta per unit cell. The results are summarized in Sec. VI. Some of the 
calculation details are presented in the Appendices \M and [B] 

II. ENERGETICALLY FAVORABLE SHAPE OF AN ISOLATED VORTEX LINE 

A. Vortex line in a finite stack of thin superconducting layers 

We start our study of the distinctive features of equilibrium vortex structures in thin films 
of anisotropic superconductors with the consideration of the vortex line shape in the layered 
systems. Let us consider a finite stack of N superconducting (SC) layers. Vortex line of an 
arbitrary shape pierces the film and can be viewed as a string of 2D pancake vortices: each 
of these pancakes is centered at the point r n = x n x + y n yo i n the n-th layer. Within the 
model of the stack of Josephson-decoupled SC layers, pancakes can interact with each other 
only via magnetic fields. We denote the interlayer spacing as s and consider each of the N 
layers as a thin film with the thickness d much less than the London penetration depth A. 
General equation for the vector potential A distribution in such system reads 



where A = X 2 /d is the effective penetration depth in a superconducting film of a vanishing 
thickness d, each n— th SC layer coincides with the plane z = z n = ns (1 < n < N), the 
sheet current at the n— th layer created by the pancake at m— th layer takes the form 



YBa 2 Cu 3 0i&&). 




(1) 



n, m=l 




(2) 
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A m (r, z) is the vector potential induced by the only pancake vortex located in the m— th 
layer (Fig. 1). The vector <fr(r) in the Eq. (j2]) is given by the expression 

$ [z x r] 



*(r) 



(3) 



2tt r 2 

and 0o = irfic/e is the flux quantum. For the layered system without Josephson coupling a 
general expression for the free energy can be written in the form: 

2 



F — — dV 
8tt 



(rotA) 2 +(^) A^J^(r)5( 



(4) 



where the total vector potential A(r, z) and the sheet current in the n— th layer J n (r), 

produced by an arbitrary vortex line are the sum of the contributions induced by all 2D 
pancakes: 

N N 

A(r, z) = J2 A m (r, z) , J n (r) = £ J™(r) . 



771=1 



m=l 



To find the magnetic vector potential A m (r, z) we adopt an approach similar to that in 



Refs. 
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equation 



22j. Between the SC layers the vector potential A m is described by the Laplace 

AA m (y)=0. (5) 



For the gauge A™ = the vector potential has only the in-plane components A m = 
(A™, A™), where 

A m (r, Z ) = -^J dqe^ r A™f/ m (q, z) , (6) 
and the function U m (q, z) can be written as 



[a™ smh. q(z n+1 — z) + a™ +1 sinhg(,2 — z n ) ] / sinh(gs) 
z n < z < z n+1 , n = 1 . . . N - 1 , 

a% exp (-q(z - z N )) , z > ztv , 



[7 m (q,^) = < 



(7) 



a™ exp — zi)) , 2; < Z\ 
Taking the Fourier transform of Eq. (j2J) we find: 



Jn(q) 



47rA 



[*(q)e* r ™<W-A™<C(q)] 



(8) 
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where 

* (q) = -j^JSM . (9) 

The sheet current density J™ results in the discontinuity of the in-plane component of the 
magnetic field B!^ across the n layer: 



Air 

— J™ = z x [Bf (r, z n + 0) - Bf (r, z n - 0)] = z 



<9A r 

Zn X 



Zn+0 

• (io) 

Zn-0 



<9;z 

Substituting the expressions ©, 0, © into above condition (flQj) we obtain the system of 
linear equations for the coefficients ot™: 

h{q) a?-a?= S lm , 

- a^_ 1 +g(q)a^-a^ +1 ^6 nm , n = 2...JV-l, (11) 

- + %) a$ = 8 Nm . 

Here we introduce two new functions which depend on the wave number q: 

g(q) = 2cosh(gs) + sinh(gs)/Ag, h(q) = cosh(gs) + (1 + 1/Ag) sinh(gs) . 

The solution of the linear system (JTT1) and the Eqs. (jSJ), (jZP define the distribution of the 
vector potential A m (r, z) which is created by a single pancake vortex positioned in the m— th 
layer. 

Without the in-plane external magnetic field Hy the relative displacement of the pancakes 
in the different layers is absent: H m k — ^m — r fe = (i.e. the pancakes form a vertical stack). 
A rather small magnetic field Hn = H a yo induces a screening Meissner current J^ 1 = J„ Xq 
in each n— th layer. Lorentz forces arising from these currents will move the pancakes from 
their initial positions. Taking into account the Eq. ( TTUj) . we find the following system of 
linear equations describing the distribution of the magnetic field screened by the layered 
structure: 



(2 + ^H 1 -H 2 = H a 



H ( > I j) 11, U,. . 0. n = 2...iV-2, (12) 



H N _ 2 +(2 + ^H N ^ = H a . 



Here H n = H n yo is the magnetic field value between the n— th and (n + 1)— th layers. The 
distribution of the Meissner screening currents in the layers takes the form: 

J?=-jf-(H 1 -H a ), J? =^-(H n -H n _ 1 ) ,n = 2...JV-2, J™ = (H a - H N ^) . 

47T 47T 47T 

(13) 

The resulting Lorentz forces F* f acting on the pancakes can be written as follows: 

Ff = (0 O / c) [Jf x z ] = (0 O / c) J„ M y . (14) 

The interaction forces between the pancakes can be found using the expression (jSJ) for the 
sheet current J™ generated by the pancake positioned in the m— th layer: 

Fr = (0o/c)[Jrxz o ] = -^-|-^5 mfc - 1 d qJl { q R mk ) a ™f 9{q) 1 g^, (15) 

where Ji(C) is the first-order Bessel function of the first kind, \ 2 ab = As = X 2 s/ d is the 
penetration depth for the in-plane currents, and 

Z(q) = 1 + 2gA/tanh(gs) . 

In order to find the equilibrium form of the vortex line in a finite stack of iV supercon- 
ducting layers under the influence of the in-plane external magnetic field Hn, we consider the 
relaxation of the set of the pancakes towards the equilibrium positions within the simplest 
version of the dynamic theory: 

V^£=J2 F n+ F n, (16) 

where rj is the viscous drag coefficient. Considering the vortex line consisting of iV = 31 
pancakes we start from the initial configuration of pancakes arranged in a straight vortex 
line parallel to the z direction (see Fig. 2). As the system approaches its final force-balanced 
(equilibrium) configuration, the velocities of all pancake motions should vanish: 

(It 

lim — - = 0, 1 < n < N . 

<->oo at 

In Fig. 2 we illustrate the evolution of the pancake configurations for two values of the 
applied in-plane magnetic field Hn and for two different numbers of layers: = 31 (Fig. 2a, 
2b) and N = 11 (Fig. 2c, 2d). The forces F^ 1 caused by the Meissner currents rotate and 

7 



bend the vortex line. For rather small applied field values this rotation and bending result 
in the formation of a certain stable configuration (see Fig. 2a, 2c). For the fields exceeding 
a certain critical value H* we do not find such stable pancake arrangement. The vortex line 
splits into two segments which move in opposite directions (see Fig. 2b, 2d). To define the 
critical value H* for the breakup of the vortex line we have carried out the calculations of 
the pancake arrangements increasing the in-plane magnetic field with the step 5H a = 0.01/Jo 
(Ho = 0o/27rA^ 6 ). The stationary vortex arrangement disappears above a certain threshold 
field value which is taken as a critical field H*. The pancake configurations for the both 
cases N = 31 and N — 11 are qualitatively similar but the values of the critical field H* for 
N = 31 and N = 11 are different. With a decrease in the number of layers N (film thickness) 
the value of the critical field H* grows: H* = 0.21H for N = 31 and H* = 0.38H for 
N = 11. In fact in layered superconductors with very weak interlayer coupling the Josephson 
vortices will appear at much lower field H a ~ <C H* . As a result, at tilted magnetic 
field, crossing lattice of pancakes, forming Abrikosov vortices, and Josephson vortices exist 
rather than a lattice of tilted vortex stacks^" 2 ^ 1 ^. The interaction between pancakes and 
in-plane field in the form of Josephson vortices produces zigzag deformation of the stack of 
the pancakes 25 . This deformation is responsible for a long range attraction between such 
stacks 2 ^ which is quite similar to the case of considered in the present work. 

B. Vortex line within anisotropic London model 

We proceed now with the consideration of the vortex line shape in an anisotropic film 
which is characterized by the London penetration depths X a b and A c for currents flowing 
parallel and perpendicular to the ab plane, respectively. We consider the case of uniaxial 
anisotropy which can be described by a dimensionless anisotropic mass tensor = mo(5ij + 
(r 2 — 1) QCj), where T = X c /X a b is the anisotropy parameter and c is the anisotropy axis. We 
choose the z axis of the coordinate system perpendicular to the film surface. In the parallel 
to the film plane direction we apply a certain magnetic field Hy = yoH a which is screened 
inside the superconducting film. 

We consider first a typical geometry when the c axis is chosen along the direction normal 
to the film plane. In such geometry the vortex line is parallel to the plane (y, z) and can 
be parameterized by a single valued function y = y(z). An appropriate thermodynamic 
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potential for determination of the energetically favorable form of the vortex line takes the 
form 

$o f D/2 , A cosh(2:/A ab ) \ 

G = F V -— dz[l , , r y {z)H a , 17 

47r J-d/2 V cosh(L)/2A afe ) J 

where D is the film thickness. The first term, F v , is the Ginzburg-Landau free energy of 

the curved vortex line, and the second term corresponds to the work of Lorentz force acting 

on the flux line and distorting this line in the presence of screening currents induced by 

the external magnetic-field component H a parallel to the film plane. To simplify the F v 

expression we consider a strong type-II superconducting material with a large ratio of the 

London penetration depths and coherence lengths. In this case the main contribution to 

the vortex line energy is determined by the energy of supercurrents j v = crotB„/47r flowing 

around the vortex core 

^ - ^/ dVrotB v (fiTotB v ) ~ ^jL- J dvfjr x ve v ,ve v ) , (is) 

where ju = m/m , and 9 V is the order parameter phase distribution around the vortex line. 
The above expression for the free energy reveals a logarithmic divergence which should be 
cut-off at both the small and large spatial length scales. The lower cut-off length is naturally 
equal to the characteristic size r c of the vortex core which is of the order of the coherence 
length. Of course, in anisotropic case one should introduce two different coherence lengths 
£ a b and £ c in the ab plane and along the c axis, respectively. The resulting core size and 
lower cut-off length r c for a certain element of the tilted flux line will, thus, depend on both 
£ a b and £ c lengths as well as on the local tilting angle of the vortex line. The upper cut-off 
length strongly depends on the ratio of the film thickness to the London penetration depth. 
For rather thick films d ^> A a & this cut-off length L r is determined by a certain combination 



of the screening lengths X a b and A c (see, e.g., Ref. 1271). For a thin film with d <C X a b one can 
separate two energy contributions: (i) the contribution coming from the region of the size 
~ d around the curved vortex line and providing the logarithmic term in the free energy 
with the upper cut-off length L c ~ d; (ii) the logarithmic contribution oc \\i(\ 2 ab /d 2 ) coming 
from the larger distances p > d which weakly depends on the vortex line shape. To sum up, 
the part of the vortex line energy which depends on its shape can be approximately written 
in the form: 



5F * ~ ik ln ^/r c ) [ m dzy/T* + y*(z) . (19) 

167r T Kb J-D/2 



Note that we neglect here the weak dependence of the logarithmic factor on the vortex line 
curvature and local orientation. Within such approximation we consider the vortex line as 
a thin elastic string which is, of course, valid provided the characteristic scale of the string 
bending is larger than the upper cut-off length L c . 

The condition of the zero first variation of the Gibbs functional gives us the following 
equation 

' w -7r*fer (20) 



u \ Ha ft cosh(z/X ab ) \ $ wr , 

b[z) = Hl t X i 1 " co 5 h (J /2A a J ■ F - = £n£ ML '^ 



where 



The equation (l20l) is valid for magnetic fields which do not exceed the critical field of the 
penetration of vortices parallel to the film plane 

I M I ^ rr(0) _ tt COsh(D/2A afc ) 

cosh( J D/2A ab ) - 1 

Thus, analogously to the case of a stack of decoupled layers the stable curved vortex lines 
can exist only for rather small magnetic fields below the critical field which corresponds 
to the penetration of a vortex parallel to the film plane. Note that in the limit H a <C H ab 



one can obtain the result found previously in Ref. 
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, TH a f cosh(z/\ ab ) \ TH a f sinh(z/A a6 ) 



H ab V cosh( J D/2A afe )7 ' tfV ' H ab V cosh(D/2A o6 ) 

Typical shape of a bent vortex line calculated from Eq. (120]) is shown in Fig. 3a. 

The above expressions can be easily generalized for an arbitrary angle x between the 
anisotropy axis c and the direction normal to the film plane. We restrict ourselves to the 
case when the axis c is parallel to the plane (y, z) and vortex line can be parameterized by 
a function y = y(z) as before. In this case the part of the free energy (jTSJ) depending on the 
vortex line shape takes the form: 

$ 2 rD/2 

n _ln(L c /r c ) / 



$ 2 r D / 2 

6F - ~ Ta 2^2 HLc/r c ) / dzy/l + y'\z) ^sin 2 0(z) + Pcos 2 0(z) 
1d7T rA^L J_ D/2 



where tan 
shape: 



(z)+x = y\ z )- Thus we find the following equation describing the vortex line 
y'(z) = + riV(,)(l + t 2 ) (22) 

y! i + t 2 r 2 (i + t 2 r 2 )v/i + t 2 r 2 -iv 2 (z) ' v ; 
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where 




t = tan x ■ 



Note that the equation (122]) is valid in the field range 



(x) = H ab cosh(D/2X ab ) / l + t 2 r 2 
cl cosh(D/2X ab )-l]l 1 + t 2 



(23) 



The critical field corresponds to the penetration of a vortex parallel to the film plane. 
Typical plots illustrating the numerical solution of the equation (122]) are shown in the Fig. 3b 
for different orientations of the applied magnetic field. Note an important difference between 
the opposite directions of the magnetic field H a : for H a > the vortex line consists of 
segments tilted in opposite directions with respect to the z axis. 

III. VORTEX VORTEX INTERACTION POTENTIAL 

In this section we derive general expressions for the interaction energy between two vor- 
tices in a thin film of an anisotropic superconductor taking into account both long range 
attraction and repulsion phenomena. We study both the limits of strong and moderate 
anisotropy for a wide range of vortex tilting angles. The shape of the interacting vortex 
lines is assumed to be fixed and not affected by the vortex-vortex interaction potential. 
Certainly, such assumption is valid only in the limit of rather larger distances between the 
vortex lines when the effect of interaction on the vortex shape can be viewed as a small 
perturbation. 

A. Interaction potential of two tilted stacks of pancakes 

In this section we consider the interaction between two vortex lines consisting of pancake 
vortices. For each vortex the pancake centers are assumed to be positioned along the straight 
line tilted at the angle 7 with respect to the anisotropy axis c (z axis). We restrict ourselves 
to the case of parallel vortex lines shifted by a certain vector R in the plane of the layers. 
Using the gauge divA = and the Fourier transform 




(24) 




(25) 
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one can rewrite the basic equation ([I]) in the momentum representation as follows: 

(q 2 + k 2 ) A(q, k) = ± £ (**(<l) ~ A M ) e ^ > ( 26 ) 



A 

n 

where 3? n (q) = $(q)e iqr ™. Taking account of the relation 



27rA„,(q) = J dke ikZn A(q, k) 



we obtain from (1261) the following equations for the Fourier components of the vector poten- 
tial A n (q): 

2gA A n = (e iqrm $(q) - A m ) e" 1 """ 11 ^ . (27) 

m 

These equations can be reduced to the scalar form 

f- + YlYl e- |n " m|9S / m = e^ r " , (28) 

m 

where we introduce the new functions / n (q): 

Jn(q) = ^ ( *„(q) - A n (q) ) = JjL *(q) / n (q) . (29) 

The solution of the linear system (1281) for a fixed distribution of pancakes r n determines 
the distribution of the vector potential A(r, z) which is created by an arbitrary vortex line 
in a finite stack of superconducting layers. In the momentum representation the general 
expression Q for the free energy of the layered system without Josephson coupling reads 

F = E / d ^ ( - A «^) ) *»(-<*) • ( 3 °) 

n J 

For two vortex lines we can write the total vector potential and the total sheet current as 
superpositions of contributions coming from the first (An , Jn ) and second (An\ J^) 
vortices: 

A n (q) = A«(q) + A^q) , J n (q) = J«(q) + j( 2 )(q). 

Calculating the interaction energy e int of vortex lines we should keep in ( 130]) only the terms 
which contain the products of fields corresponding to different vortex lines: 

= E / rf M ( *? } (q) " A«(q) ) < 2) (-q) + ( $i 2) (q) " A^(q) ) *«(-q) ] . 

(31) 
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Finally, for the particular case of two parallel vortex lines which are shifted at the vector 
R = r„ — rj^ (n = 1, N) in the (xy) plane we get following expression for the interaction 
energy via the scalar functions / n (q) : 

£ int(R) = Sx cos(qR) £ /»(q) • (32) 

The expression (132]) and equations (j28|) determine the interaction energy of two identically 
bent vortex lines. 

Our further consideration in this subsection is based on two assumptions: (i) for each 
vortex we choose the centers of pancakes to be positioned along the straight line tilted at a 
certain angle 7 relative to z axis, and put r n = rastan7yo; (ii) we restrict ourselves by the 
continuous limit assuming qs and q y s tan 7 <C 1. In this case the Eqs. (1281) . (|32|) can be 
simplified (see Appendix lAl for details): 

e mt (R) = —4— / d\ cos(qR) S(q) , (33) 
S{q) = t ~ \ D- 



Aab q 2 { l + p 2 

^2(1 — k 2 ) [k(l — p 2 ) sinhL + (k 2 — p 2 )(coshL — cos(pL)) + 2kpsin(pL)] . t ^ ^ , 
U 2 + Kb (1 + P 2 ) 2 t 2fc c osh L + (1 + A; 2 ) sinh L)] 



where 



L = D^q* + \- 2 , k = q/ y q 2 + \~ 2 , p = q y tan 7/^ + Kb , (35) 

and D = (N — l)s is the thickness of the superconducting film. The first term in (|34|) 
describes the interaction in the bulk system, while the second term is responsible for the 
effect of film boundaries. 

The minimum energy configuration corresponds to the case R x = 0. In Fig. 4 we present 
some typical plots of the interaction energy e int (R x = 0,R y ) vs the distance R y = R for 
d = 3X a b which corresponds to the Lorentz microscopy experiments in YBCO^ and Bi-2212^ 
samples. Analyzing the dependence Ei n t(R), one can separate three contributions to the 
energy of vortex-vortex interaction: (i) a short-range repulsion which decays exponentially 
with increasing intervortex distance R (for R > \ a b); (ii) an intervortex attraction which 
is known to be specific for tilted vortices in anisotropic systems; this attraction energy 
term decays as R~ 2 and strongly depends on the angle 7 between the vortex axis and the 
c direction; (iii) long-range (Pearl) repulsion which decays as R' 1 and results from the 
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surface contribution to the energy. Note that the third term does exist even for a large 
sample thickness D (see Ref.— ) although in the limit D ^> X a b it is certainly masked by the 
dominant bulk contribution. At R ^> A & the short-range interaction term vanishes and the 
interaction energy vs R takes the form 

^ 4 f D e// tan 2 7 2 
m * " 8tt 2 V # 2 R 

where D e ff = D — 2X a b tanh(D/2\ a b) is the effective film thickness. One can observe here an 
interplay between the long-range attractive (first term in Eq. fl36p ) and the repulsive (second 
term in Eq. (l36]) ) forces. Note that the X a b value increases with an increase in temperature, 
thus, the effective thickness decreases and the long range attraction force appears to be 
suppressed with increasing temperature. For large R the energy is always positive and 
corresponds to the vortex repulsion similar to the one between the pancakes in a single layer 
system. With a decrease in the distance R the attraction force comes into play resulting in 
the change of the sign of the energy. Such behavior points to the appearance of a minimum 
in the interaction potential. 




B. Vortex vortex interaction within anisotropic London model 

We now proceed with the derivation of the expression for the intervortex interaction en- 
ergy in an anisotropic superconducting film. We choose the anisotropy axis c (z— axis) to be 
oriented perpendicular to the film plane and consider two curved vortex lines with identical 



shapes 
in Ref. 



bund in Sec. II B. Our further calculations are based on general expressions derived 
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for the energy of an arbitrary vorticity distribution in an anisotropic supercon- 
ducting film of finite thickness (see Appendix[B]for details). The resulting interaction energy 
of two curved vortices shifted from each other in the y direction at a certain distance R can 
be presented in the form: 

£int = £0 {^nt + ^nt^ + e int) > (37) 

where eo = 0o/167r 3 A a 6, while e v int , el n r t ay , and e^ c are given by the expressions ( 1B5I) . ( 1B6I) . 

(EzD. 

Considering the particular case of straight vortex lines parallel to the plane (yz) and 
tilted at a certain angle 7 with respect to the c direction we obtain the following expression 
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for the interaction energy of two vortices: 

e int (R) = — [ d\ cos(qR) S r (q) , (38) 
1 f „ fl + A# 1 



Sr(q) = t— j ^ 



Aatg 2 I Vi + Pr (i + g 2 )(i+P 2 ) 

2(1 — A; 2 ) [£;(1 — p 2 ) sinhL + (fc 2 — p 2 )(coshL — cos(pL)) + 2kpsin(pL)} 
q 2 + A" 2 (1 + p 2 ) 2 [2k cosh L + (1 + A; 2 ) sinh L)] 



,(39) 



where 



fc r = g tan 7/ \/T 2 g 2 + X ab , p T = q y tan 7/ y r V + A afe , 

and the parameters L, k and p are described by the expressions (|35|) . In the limit of 
strong anisotropy (r >> 1) the spectral function S*r(q) fl39l) naturally coincides with the 
corresponding expressions (13"1|) obtained for the layered system without Josephson coupling. 

Some typical plots of the interaction energy vs the intervortex distance for tilted vortex 
lines calculated using the Eqs. fl38|) .f l39|) are shown in Figs. I5|6l Analyzing the dependence 
£int{R) one can separate three contributions to the energy of intervortex interaction: (i) a 



short-range repulsion (for R A^a/I + tan 2 7) which decays exponentially with increasing 
intervortex distance R; (ii) an intervortex attraction which comes into play for the region 



AabV^l + tan 2 7 < R < T\ a b and decays exponentially with the vortex-vortex distance R 
for R > TXaj,; (iii) long-range (Pearl) repulsion which decays as i? -1 at large distances and 



results from the surface contribution to the energy. Taking the limit R A^^l + tan 2 7 
we get 



. D7r v /r 2 + tan 2 7 1 (L c 
£int/£o ^ ^ \~R 



In the region A^yl + tan 2 7 < R < TX^ the short-range interaction term vanishes and 
the interaction energy vs R is given by the sum (!36l) of attractive and Pearl's contributions. 
Similarly to the case of decoupled layers discussed above the attractive term can result in the 
appearance of a minimum in the dependence of the vortex-vortex interaction potential vs R. 
The position of this minimum can be roughly estimated as the boundary of the region of the 



short-range repulsion: R m i n ~ A^a/I + tan 2 7. Obviously, the minimum should disappear 
provided R m in > TAab, i-e., when the region of the attraction between vortices vanishes. 
This condition gives us the the upper boundary on the tilting angle 7 restricting the interval 
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of the energy minimum existence: 

tan 2 7 < T 2 - 1 . 

The lower boundary of this angular interval can be found comparing the attractive and 
repulsive terms in the expression fl36l) at the distance R m in- 



These analytical estimates of the angular interval are in a rough qualitative agreement with 
the numerical calculations (see Figs. EE]) for two values of the film thickness D = 3X a b, 10A a b. 
Indeed, one can see that increasing the tilting angle we first deepen the minimum in the 
interaction potential and then make it more shallow. The figures confirms the deepening of 
the minimum with the increase in the anisotropy parameter T. Our numerical calculations 
demonstrate that for the film thickness D = 3X a b (Fig. [5]) the minimum of the interaction 
energy of two straight tilted vortices can appear only for T > 14. Starting from r w 27 
the bound vortex pair becomes energetically favorable. An increase in the film thickness 
reduces the relative contribution of Pearl repulsion to the energy of intervortex interaction 
Bint- As a result attraction of vortices takes place for smaller values of the tilting angle 7 and 
anisotropy parameter V. Thus, in a film with the thickness D = 10A a {, (Fig. [6]) the minimum 
in the ei nt (R) dependence appears for T > 7, whereas creation of the bound vortex pair 
becomes energetically favorable for T > 9. 

As a next step, we check if the above results obtained for straight tilted vortices remain 
valid for the curved vortex lines. Our analysis of the effect of the vortex line curvature is 
carried out for model vortex profiles found in Sec. II B. The resulting typical dependencies of 
the intervortex interaction potential vs R for different magnetic field values and anisotropy 
parameters are shown in Figs. 171(81 One can clearly see that the minimum in the interaction 
potential vs R survives when we take account of the vortex line curvature. Moreover the 
curving of the vortex line even deepens this minimum as it is confirmed by the comparison 
of energies of straight tilted and curved vortices presented in Fig. [HJ For such comparison 
we choose the straight vortex lines connecting the ends of curved vortices. We find that 
for curved vortices the energy minimum exists even for smaller anisotropy parameters than 
for straight vortices (i.e., the threshold anisotropy value for D = 3X a b becomes less than 




16 



r ~ 14). Of course, increasing the film thickness one can weaken the restrictions on the 
existence of the minimum in the interaction potential: e.g., for D = 10A a & the minimum 
appears at T > 9. 

The above theoretical analysis demonstrates that vortex-vortex attraction and the for- 
mation of chains are possible only for the rather large tilting angles and at low vortex 
concentrations, i.e., when the magnetic-field component H z perpendicular to the film plane 
is very weak. In fields H z slightly above H c \ Abrikosov vortices will form chains due to 
the long range attractive interaction. Peculiarities of penetration of such chains of tilted 
Abrikosov vortices into bulk layered (anisotropic) superconductor are well known: in the 
first approximation, the vortex period in chains does not depend on applied magnetic field, 
while the distance between chains changes as 1/H Z . The presence of vortex chains sig- 
nificantly modifies the magnetization curves with respect to analogous curves for isotropic 
superconductors.— In the next sections we discuss additional peculiarities of intervortex 
interaction specific for thin-film samples of layered (anisotropic) superconductors. 

IV. VORTEX CLUSTERS 

The unusual vortex-vortex interaction potential behavior discussed in the previous sec- 
tion can result in unconventional vortex structures. We start our analysis of energetically 
favorable vortex structures from the problem of stability of a vortex chain. The formation 
of infinite vortex chains is known to be a signature of the intervortex attraction in bulk 
anisotropic superconductors. The long range repulsion of vortices in thin films can destroy 
the infinite vortex chains. Indeed, despite of the fact that two vortices attract each other 
at a certain distance, further increase in the number of vortices arranged in a chain can be 
energetically unfavorable because of the slower decay of the repulsive force compared to the 
attractive one. In this case, for rather thin samples, there appears an intriguing possibility 
to observe vortex chains of finite length, i.e., vortex molecules or clusters. In this section we 
present the calculations of energies of such vortex clusters. 

As we have demonstrated above, the minimum in the interaction potential exists for both 
the limits of strong and moderate anisotropy. The vortex molecule cohesion energy is given 
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by the expression: 

^ = J2B int (R ij ) (40) 

i>j 

where N is the number of vortices in the molecule, and Rij are the distances between 
i— th and j—th vortices in the chain molecule. Shown in Figs. l9][T0l are typical plots of the 
interaction energy per vortex vs the intervortex distance R for equidistant vortex chains 
with different N numbers calculated within the model of decoupled superconducting layers 
and anisotropic London theory. The energetically favorable number of vortices in a molecule 
grows as we increase the film thickness and/or the tilting angle because of the increasing 
attraction term in the pair potential 8^. Shown in the insets of Figs. [9]are schematic pictures 
of vortex matter consisting of dimeric and trimeric molecules. Finally, for rather thick 
samples with D ^> A a & we get a standard infinite chain structure typical for bulk systems. 
Note that the formation of an infinite vortex chain may be considered in some sense as a 
polymerization of the vortex molecules. Certainly, the crossover from the vortex molecule 
state to the infinite chain structure is strongly influenced by the increase in the vortex 
concentration governed by the component B z of the external magnetic field perpendicular 
to the film. Indeed, one can expect such a cross-over to occur when the mean intervortex 
spacing approaches the molecule size. Thus, the vortex molecule state can appear only in a 
rather weak perpendicular field when its observation can be complicated, of course, by the 
pinning effects. 

V. PHASE TRANSITIONS IN VORTEX LATTICES 

Considering the effect of a finite magnetic field (i.e., a finite concentration of vortex 
clusters) we restrict ourselves by the simplest case of regular vortex arrays. For a regular 
vortex array the formation of clusters corresponds to the transition with a change in the 
number of vortices in the elementary lattice cell. The mechanism underlying such transition 
is naturally connected with the appearance of the minimum in the interaction potential 
for a vortex pair. In this section we present our calculations of energy of vortex lattices 
with different number of flux quanta per unit cell. The possibility to get the energetically 
favorable states with a few vortices per unit cell will be illustrated for a particular intervortex 
interaction potential derived above for a model of decoupled superconducting layers. The 
generalization of such consideration for anisotropic London theory is straightforward. Note 
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that the vortex lattice structure for bulk anisotropic superconductors in tilted field in the 
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framework of London approach has been calculated in Ref. 

Let's consider a vortex lattice characterized by the translation vectors T = zai + ja 2 , 
where a 12 are primitive vectors of the lattice. The primitive cell occupies the area A Q = 
[&i x a 2 ] • z and is assumed to contain M vortices: B Z A = M$ - Positions of vortices in a 
cell are determined by the vectors r m (m = 1, M) (see Fig. [TT]) . The interaction energy per 
unit lattice cell can be expressed via the vortex-vortex interaction potentials fl33|) . fl34|) : 

M M 
^c( r mfc,T)= £int( r mk) + Ej nt (T + r m k) ■ (41) 

m, ky^m T^O m, k 

The interaction energy ( I4ip depends on both the relative positions r m fc = r m — of vortices 
in the primitive cell and the structure of the vortex lattice defined by the translation vectors 
T. The first term in ( 14T1) describes the interaction energy between vortices in the primitive 
cell (without the lattice contribution), whereas the second sum takes account of the lattice 
effects. With the help of the Poisson formula, one can rewrite the intervortex interaction 
energy ( T4T]) in terms of the Fourier components 



16vr 3 A 



06 



A 9 M 

47T 2 



,2 <" r 

-^r^TS(Q) cos(Qr mfc )-M J d 2 qS(q) 



A, 

Q m, k 



(42) 



where the function S(q) is determined by the Eq. fl34|) . and Q are the reciprocal-lattice 
vectors. The sum and the integral in Eq. (14"2"]1 diverge both at Q = and at large Q values. 
The small Q divergence corresponds to the linear (in the system size) increase in the vortex 
energy because of the slow 1/R decay of the vortex- vortex interaction potential. The large 
Q divergence is logarithmic and is associated with the vortex self energy. 

For simplicity, we restrict ourselves to the case of an instability with respect to the unit 
cell doubling and tripling, i.e. formation of the vortex lattices with two and three flux quanta 
per unit cell (M = 2 and M — 3). Hereafter we consider only the shifts of vortex sublattices 
along the y direction and choose the appropriate reciprocal-lattice vectors 

27T 71 

Qij = -j- (* - J/4) x + - j y , = 0, ±1, ±2, . . . , 

n 2tt . . , 2vr 

Qij = y (* - J/ 6 ) x o + 3- J yo , 1 13 — 0, ±1, ±2, . . . 

for M = 2 and M = 3, respectively. Here we consider only equidistant vortex chains within 
the primitive cells. Fixing the value of the field B z we fix the unit cell area area A = 2ab 
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for M = 2 and Aq = 3ab for M = 3. Thus, the interaction energy f|4"2]) depends only on two 
parameters: (i) a = b/a ratio characterizing the lattice deformation; (ii) relative displace- 
ment Aa of vortex sublattices along the y-axis (see Fig. [TTjl . To exclude the divergence at 
Q = it is convenient to deal with the energy difference: 



The results of our numerical calculations of this energy difference are shown in Fig. [12k . 
One can clearly observe that changing the vortex tilting angle we obtain the minimum in 
the function Ae c (Aa) which gives us the evidence for the phase transition in the lattice 
structure with the unit cell doubling or tripling depending on the vortex tilting angle. The 
multiplication of the unit cell is accompanied by the strong change in the lattice deformation 
ratio a (see Fig. [T2"b). 

VI. CONCLUSIONS 

To sum up, we suggest a scenario of the phase transitions between the flux structures with 
different number of vortices per unit cell which can be realized in thin films of anisotropic 
superconductors placed in tilted magnetic fields. We demonstrate that the vortex inter- 
action in the films of anisotropic superconductors placed in tilted magnetic fields is very 
special. The underlying physics arises from the interplay between the long range attrac- 
tion and repulsion between tilted vortex lines. In consequence, new and very reach types 
of vortex structures may appear. They are formed from the vortex dimers, trimers, etc., 
and the transition between different types of vortex structures may be controlled by tilting 
of external magnetic field and/or by varying of the temperature. Our theoretical findings 
are based on two theoretical approaches: anisotropic London model and the London-type 
model of decoupled superconducting layers. Taking account of the vortex tilt and bending 
we analyzed the distinctive features of the vortex-vortex interaction potential in a wide 
range of parameters and fields and demonstrated the possibility to obtain a minimum in 
the vortex interaction potential vs the intervortex distance. Further analysis in the paper 
included the calculations of energy of finite size vortex clusters as well as the energy of 
regular vortex arrays with different number of vortices per unit cell. The phase transitions 
accompanied by the multiplication of the primitive lattice cell appear to be possible for 




(43) 



(J 



<7 
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dilute vortex arrays, i.e. for rather small magnetic field component B z . We believe that 
our theoretical predictions concerning the unusual vortex configurations are experimentally 
observable using the modern vortex imaging methods such as Lorentz microscopy, scanning 
tunneling microscopy, scanning Hall-probe or decoration technique. 
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Appendix A: Interaction energy of tilted vortices: continuous limit 

Let us evaluate the interaction energy (1321) of two tilted parallel vortex lines taking 

= ns tan 7 x , r {2) = r^ 1} + R 

and assuming qs <C 1 and q x stwa. r j <C 1. We introduce a continuous coordinate t = ns and 
continuous function / q (£). Thus, the linear system of equations (|2"5|) reduces to the following 
integral equation 

D/2 

AW + tV / dt'e^-^Ut 1 ) = e^' ta ^ . (Al) 

l( l X ab J 

-D/2 

The equation (lAlj) can be rewritten as a differential one 

% ~ (Kb + <?) U(t) = - (<& tan 2 7 + q 2 ) (A2) 
at the interval —D/2 < t < D/2 with the boundary conditions 

'dfq 



dt 

Introducing the notations 



±9/q 



= (ifetan 7 ± _ ^ 

±D/2 



T = tJq 2 + \- 2 , L = DJq* + \- 2 , k = q/Jq 2 + \7 2 p = q x tan 7/ Jq 2 + A 



-2 
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one can rewrite the equation (1A2I) and boundary conditions (1A3I) in dimensionless form 

*% - /, = - (P 2 + k 2 ) , (A4) 



= (ip ± A;)e ipL/2 . (A5) 

±L/2 



The solution of the Eq. (IA4[) has the form 



p 2 + /c 2 
i + P 2 

where the constants a and b are defined by the boundary conditions flA5j) : 

_ e ( ip+1 ) L / 2 (fc + + fc) + e"^ +1 ) L / 2 (fc - ip)(l - k) 
a ~ 2(2£; cosh L + (1 + k 2 ) sinh L) 

_ e (ip-l)V2(jfc + _ jfc) + e (-*+i)x/2^ _ ip ^! + fc ) 

~ 2(2/c cosh Z + (1 + A; 2 ) sinh Z) ' 

In the continuous limit the expression for the interaction energy (1321) takes the form: 

*3 /" ^ 2 q 



Here the function 

D/2 



^ = 77~v\ir / -r cos(qR) 5(q) . (A7) 



5(q)= y rft/ q (t)e- ifc * tan7 



-D/2 

can be calculated analytically: 

_ . _ ^p 2 + fc 2 ^2(l - k 2 ) [k(l - p 2 ) sinhL + (k 2 -p 2 )(coshL - cos(pL)) + 2A;psin(pL)] 
1+p2 ^/g 2 + A a6 2 (l+p 2 ) 2 [ 2k cosh £ + (l + jfe 2 ) sinh L] 

(A8) 

Appendix B: Interaction energy of curved vortices: anisotropic London model 

To calculate the vortex-vortex interaction within anisotropic London model we use gen- 



eral expressions derived in Ref. 
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for the total energy E of an arbitrary arrangement of 
curved vortices in a superconducting film of thickness D with the c-axis perpendicular to 
the film plane: 

E = E vi + E stray + E vac , (Bl) 

where 
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2 

my I j 



E 



stray 



1 /■ d 2 k i A; 2 



87T ./ 47T 2 V 



(l-e- 2 ° D )h+| 2 + (e 2 » D -l)| 7 -| 2 ) 



^=s/^( e ^ D l*-| , + l* + l') 



Here a = x,y,z and 
G x (k) = G„(k) = 



, G,(k) 



1 + k 2 X 2 r 



1 + A; 2 A? + k 2 \ 2 ab ' ZW (1 + ^(1 + k 2 ± X 2 + k 2 \ 2 J ' 



v 



7"(k±) 



fciC 



, 7 + (kx) 



fc,C 



r(k ± ) = (v/A^C) x 
0+(k ± ) = (-v/k ± C) 



2k ± B + 



2k i A 



(k ± + v)e vD + (k ± -v)e 



(k ± + v)e vU + (k ± - v)e 



-vD 



-vD 



A 



B 



C(k ± ) = e~ vD (k ± - v) 2 - e vD (k ± + vf , v = Jk\ + \J , 



-2 , \~2 



9i(k) 



4>v, 



1 + k 2 \ 2 ab ■ 



The summation in the above expressions is carried out over k z = k m = mn/D, m = 
0, ±1,±2, ... and _L stands for the vector component parallel to the xy plane. Following 
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Ref. |29j we introduce here the Fourier transform u m (k±, k m ) of the vorticity distribution 



u"'[T): 



v± m (k±,k m ) = -2% \ d 2 r L e 



u z vi (k ± ,k m ) = 2 / d 2 r ± e 



-ik i r i 



-ik j_ rj 



dz 8m(k m z)v™(vi_, z) 



D 

o 

D 



dz cos(k m z)i>™(r±, z) . 



(B2) 
(B3) 



For a pair of curved vortices shifted in the y direction at a certain distance R the expressions 
(El and (El take the form: 



-2i I 1 + e~ ikyR 



is™{k x ,k m ) = 



D 



v v z l (k ± ,k m ) = 2 l + e 



D 



To find the vortex-vortex interaction energy we should take the terms in Eq. (IBlj) which 
depend on mutual vortex arrangement: 



_ / i/i , stray . V ac\ 
tint — fc V, fc int ' & int fc mt J 



where e — 0o/16vr 3 A a 6 and 



(B4) 



poo pO pO f 

Ct = *] vdvj dd J d( 2 I Jo [v (pid) - tKCi) + 

x(n 1 ( I /,Ci,C 2 ) + n 2 (z/,Ci,C 2 )) , 



^stray 
-int 



4tt / vdv \ d(i / dC2 < Jo 



D J-D 



Jn 



^h?(Ca)-*KCi)+£ 



+ Jo 



X- 



n 3 (^Ci,C 2 ) 



(B5) 

" [viCi) - v(Ci) - R 

(B6) 



r(z/) sinh r(z/) J D <^ e - T M 6 [u - r(z/)] 2 - e T ^ 3 [u + r(u)f 



2 • 



vac 



dCi d( 2 < Jo 

n 4 (^,Ci,C 2 ) 



v(ti(C2)-v(Ci) + R) +Jo v(v(<Z2)-ri(Ci)-R 



z/sinh 



t(u)D I e-<")» [v - r(u)) 2 - [v + t(u)\ 



2 ! 



(B7) 
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r(u) = Vl + TV 2 , u{u) = Vl + v 2 . 

Here we have introduced dimensionless coordinates rj = y/\ a b, C — z /^ab, dimensionless 
wave number v = q\ a b, and use the notations 



n 1 (z/,Ci,C 2 )=V(Ciy(C 2 ) 



sinh 


t(iz) (£> + (Ci + C2 - ICi - C2D/2) 


" 8inh[r(i/)((-Ci-C2-|Ci-C2|)/2)] 


t(v) sinh 







n 2 (z/,Ci,C 2 ) = 



t{v) cosh 


t(i/) (£> + (Ci + C 2 -|Ci-C2|)/2) 


' oosh[r(i/)((-Ci-C a -|Ci-C2|)/2)] 


cosh 


z/ 2 sinh 

u(u) (d + (Ci + C2-|Ci-C2|)/2)' 


Z>t(i/) 
cosh [w 


M((-Ci-C 2 -|Ci-C2|)/2)] 



v 2 u[y) sinh Du(v) 



n 3 (i/,Ci,C2) = <^ 2 cosh r(i/)D +r(z/) 2 cosh r(i/)D + 2i/r(i/) sinh r(i/)D 



r(i/)(D + 0) cosh r(D + C 2 ) + cosh [r(z/)Ci] cosh [r(z/)0] 



x < cosh 



+ cosh 



t(is){D + Ci)]) cosh [r(z/)C 2 ] + cosh [t(u) Zi ] cosh r(i/)(C 2 + £>) 



114(^,(1,(2) = i 2i/ 2 cosh 2 r(i/)D + sinh 2 t(u)D + 2z/t(z/) sinh t(z/)L> cosh r{u)D 



x 



cosh 



r(i/)(D + Ci) cosh t(v)(D + ( 2 ) + cosh [r(z/)Ci] cosh [r(z/)d] 



2<u cosh 



t(i/)L> 



+ vt sinh 



r(i/)£>j [> x 



cosh T{y){p + Ci)J cosh [r(z/)C 2 ] + cosh [r(i/)Ci] cosh t{v)(C, 2 + L>) 

The dimensionless thickness of the film D and dimensionless intervortex distance R are 
measured in the units of A a &. 
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FIG. 1: A single 2D pancake vortex positioned in the m— th layer of a finite layered structure, d is 
a thickness of the superconducting layer, and s is the distance between the layers. 
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FIG. 2: (Color online) Configurations of N = 31 (panels a, b) and iV = 11 (panels c, d) pancakes 
in a finite stack in the presence of the applied in-plane magnetic field H a . (a) The force-balanced 
(equilibrium) configuration of pancakes for H a = 0.2Hq < H* . (b) Pancake configurations at 
sequential time points t\ < t<i < t% for H a = 0.22//o > H* . For the structure with N = 31 we find 
H* ~ 0.21i?o- (c) The force-balanced (equilibrium) configuration of pancakes for H a = 0.35i7o < 
H*. (d) Pancake configurations at sequential time points ii < ti < t% for H a = OAHq > H*. For 
the structure with N = 11 we find H* ~ 0.38i? - Here H = o /2ttA^, A = 10A afe , and s = 0.1A o6 . 
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FIG. 3: (Color online) Typical configurations of the vortex lines in the film of the thickness 
D = 3A a t for the anisotropy parameter T = 5 and for different values of in-plane magnetic field 
H\\ = H a yQ. (a) The anisotropy axis is perpendicular to the film plane (x = 0°, ~ 1.74 H a f,). 
(b) The anisotropy axis is tilted with respect to the z axis (x = 30°, Hj$ — A.QH^) The numbers 
near the curves denote the values of the ratio H a /H a i,. The dashed line shows the shape of a vortex 
line in the absence of the in-plane magnetic field. 




FIG. 4: (Color online) Typical plots of the interaction energy per vortex [Eqs. (I33p and (!34p ] vs 
the distance R between two vortices for a film of thickness d = 3\ a b and different tilting angles 
7 = 70°, 75°, 78°, 80° (e = ^/16vr 3 A ab ) . 
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FIG. 5: (Color online) Typical plots of the interaction energy per vortex [Eqs. ([38j) . ([39jl ] vs the 
distance R between two tilted vortices for an anisotropic film of the thickness D = 3\ a b- (a) 
Interaction energy for the anisotropy parameter T = 27 and different tilting angles. The numbers 
near the curves denote the values of tilting angle 7. (b) Interaction energy for 7 = 83° and different 
values of anisotropy parameter. The numbers near the curves denote the values of T. 




FIG. 6: (Color online) Typical plots of the interaction energy per vortex [Eqs. ([38j) . (f3"9]) ] vs the 
distance R between two tilted vortices for an anisotropic film of the thickness D = 10A a {,. (a) 
Interaction energy for the anisotropy parameter T = 9 and different tilting angles. The numbers 
near the curves denote the values of tilting angle 7. (b) Interaction energy for 7 = 73° and different 
values of anisotropy parameter. The numbers near the curves denote the values of T. 
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FIG. 7: (Color online) Typical plots of the interaction energy per vortex (|37p vs the distance R 
between two curved vortices for an anisotropic film of the thickness D = 3X a b'. (a) T = 15 ; (b) 
r = 27. The numbers near the curves denote the values of the ratio H a /H a b. The shape of vortex 
lines is schematically shown in the insets. 




FIG. 8: (Color online) Comparison of the vortex-vortex interaction potentials for curved [Eq. (I37p ] 
(solid lines) and straight tilted [Eqs. ()38|) . (|39p ] (dashed lines) vortices for an anisotropic film of the 
thickness D = ?>\ a b with different anisotropy parameters: (a) V = 15, H a = 0.91i^ a t (7 = 80.6°); 
(b) r = 27, H a = 0.6H a i, (7 = 81.5°). The shape of vortex lines is schematically shown in the 
insets. 
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FIG. 9: (Color online) Typical plots of the interaction energy per vortex [Eqs. (f33}) . (fMj) . (fi0j) ] vs the 
intervortex distance R in an equidistant chain of N vortices in a stack of decoupled superconducting 
layers (D = 3\ a b): (a) 7 = 78°; (b) 7 = 80°. The numbers near the curves denote the number N 
of vortices in molecule. Inserts show schematic pictures of vortex matter consisting of dimeric (a) 
and trimeric (b) molecules. 
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FIG. 10: (Color online) Typical plots of the interaction energy per vortex vs the intervortex 
distance R in an equidistant chain of N vortices for d = 3X a b, T = 27, H a = 0.513-£f a b(7 = 80°). 
The numbers near the curves denote the number N of vortices in a molecule. The shape of vortex 
line and effective tilting angle 7 are schematically shown in the inset. 
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FIG. 11: (Color online) Vortex lattice with two M = 2 (a) and three M = 3 (6) vortices per a 
primitive cell. 




FIG. 12: (Color online) (a) The energy difference Ae c vs the relative displacement Aa of vortex 
sublattices for different tilting angles 7 = 78° (solid line) and 7 = 80° (dashed line) and different 
number of flux quanta per unit cell M = 2, 3. (b) Lattice deformation ratio a = b/a vs the relative 
displacement Aa of vortex sublattices for different tilting angles 7 = 78° (solid line) and 7 = 80° 
(dashed line) and different number of flux quanta per unit cell M = 2, 3. Here we put ao = 60A a b. 
The numbers near the curves denote the number M of vortices per unit cell. 
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